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Abstract A configuration p in r-dimensional Euclidean space is a finite collection 
of labeled points p , . . . ,p n in W that affinely span W. Each configuration p defines 
a Euclidean distance matrix D p = (d,j) = (\\p l — p-*\\ 2 ), where || || denotes the Eu- 
clidean norm. A fundamental problem in distance geometry is to find out whether 
or not, a given proper subset of the entries of D p suffices to uniquely determine the 
entire matrix D p . This problem is known as the universal rigidity problem of bar 
frameworks. In this chapter, we present a unified approach for the universal rigidity 
of bar frameworks, based on Euclidean distance matrices (EDMs), or equivalently, 
on projected Gram matrices. This approach makes the universal rigidity problem 
amenable to semidefinite programming methodology. Using this approach, we sur- 
vey some recently obtained results and their proofs, emphasizing the case where the 
points p 1 , . . . ,p n are in general position. 



1 Introduction 

A configuration p in r-dimensional Euclidean space is a finite collection of la- 
beled points p l ,...,p n in W that affinely span W. Each configuration p defines 
the nxn matrix D p = (dij) = (\\p l — p-*\ | 2 ), where 1 1.| | denotes the Euclidean norm. 
D p is called the Euclidean distance matrix (EDM) generated by configuration p. 
Obviously, D p is a real symmetric matrix whose diagonal entries are all zeros. A 
fundamental problem in distance geometry is to find out whether or not, a given 
proper subset of the entries of D p , the EDM generated by configuration p, suffices 
to uniquely determine the entire matrix D p ; i.e., to uniquely recover p, up to a rigid 
motion. This problem is known as the universal rigidity problem of bar frameworks. 

A bar framework, or framework for short, denoted by G(p), in W is a configu- 
ration p in W together with a simple graph G on the vertices 1,2, ... ,n. To avoid 
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trivialities, we assume throughout this chapter that graph G is connected and not 
complete. It is useful to think of each node i of G in a framework G(p) as a univer- 
sal joint located at p\ and of each edge (2, j) of G as a stiff bar of length | \p l — p-> 1 1 . 
Hence, a bar framework is often defined as a collection of stiff bars joined at their 
ends by universal joints. Figure Q] depicts a framework G(p) on 4 vertices in K 2 , 
where G is the complete graph K4 minus an edge, and the points p l , . . . ,/? 4 are the 
vertices of the unit square. 

4 3 
a ) a 

y- ^ 2 

Fig. 1 A bar framework G(p) on 4 vertices in R 2 , where V(G) = {1,2,3,4}, E(G) = 
{(1,2), (2,3), (3,4), (4, 1), (1,3)} and p 1 ,p 2 ,p 3 ,p 4 are the vertices of the unit square. 

We say that two frameworks G(p) and G(q) in W are congruent if D p = D q . 
Furthermore, let H denote the adjacency matrix of graph G, then two frameworks 
G(p) in W and G(q) in W are said to be equivalent if H oD p = H oD q , where o de- 
notes the Hadamard product, i.e., the element- wise product. We say that framework 
G(q) in W is affinely-equivalent to framework G(p) in W if G(q) is equivalent to 
G(p) and configuration q is obtained from configuration p by an affine motion; i.e., 
q' =Ap' + b, for all 1 = 1 , . . . , n, for some r x r matrix A and an r-vector b. 

A framework G(p) in W is said to be universally rigid if every framework G(q) 
in any dimension that is equivalent to G(p), is in fact congruent to G(p); i.e., if for 
every framework G(q) in any dimension such that H o D q = H o D p , it follows that 

Dq=D p . 

Thus, given D p = (djj), the EDM generated by configuration p, let C '■ 
i < /'; for ij — 1,2, ...,«}. Then the proper subset of entries of D p given by 
: (z',y) G /T} suffices to uniquely determine the entire matrix D p if and only if 
framework G(p) is universally rigid, where G — (V,E) is the graph with vertex set 
V = {1 ,2, . . . ,re} and edge set E =K. For example, the framework given in Figure[T] 
is not universally rigid; and the subset of entries of D p given by {dij : (i,j) £ £(G)} 
does not uniquely determine the entire matrix D p since the entry ^24 can assume 
any value between and 2. 

The notion of dimensional rigidity is closely related to that of universal rigidity. 
A framework G(p) in W is said to be dimensionally rigid if there does not exist 
a framework G(q) that is equivalent to G(p), in any Euclidean space of dimension 
> r + l. For example, the framework G(p) given in Figure Q] is obviously not di- 
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mensionally rigid since there is an infinite number of frameworks G(q) in M? that 
are equivalent to G(p). 

In this chapter, we survey some recently obtained results concerning framework 
universal as well as dimensional rigidity. These results are given in Section [2] and 
their proofs are given in Section|4] Section[3]is dedicated to the mathematical pre- 
liminaries needed for our proofs. Our EDM approach of universal rigidity of bar 
frameworks extends to the closely related notion of "local" rigidity. However, due 
to space limitation, "local" rigidity [3 | will not be considered here. Also, we will 
not consider the other closely related notion of global rigidity lfT0l[T3ll . 



2 Main Results 

The following theorem characterizes universal rigidity in terms of dimensional 
rigidity and affine-equivalence. 

Theorem 1 (Alfakih |2|). Let G(p) be a bar framework on n vertices in W, r < 
n — 2. Then G(p) is universally rigid if and only if the following two conditions 
hold: 

1. G(p) is dimensionally rigid. 

2. There does not exist a bar framework G(q) in W that is affinely-equivalent, but 
not congruent, to G{p). 

The proof of Theorem[T]is given in Section|4] The notion of a stress matrix S of a 
framework G(p) plays an important role in the characterization of universal rigidity 
of G(p). Let G(p) be a framework on n vertices in R'', r < n — 2. An equilibrium 
stress of G(p) is a real valued function co on E(G), the set of edges of G, such that 



£ (Oijitf - p j ) = for all i= 1,. 

j:(iJ)GE(G) 



.,n. 



(1) 



Let co be an equilibrium stress of G(p). Then the nxn symmetric matrix S = 
where 

-coij if (U) eE(G), 

if*Vjand(U)££(G), 
£ (Oik if i = j, 

t.(i,k)€E(G) 



(2) 



is called the stress matrix associated with co, or a stress matrix of G(p). 

Given framework G(p) on n vertices in W, we define the following nx r matrix 



,2 r 



P:= 
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P is called the configuration matrix of G(p). Note that P has full column rank since 
p , . . . ,p" affinely span W. The following lemma provides an upper bound on the 
rank of a stress matrix S. 

Lemma 1. Let G(p) be a bar framework on n nodes in W, r < n — 2, and let S and 
P be a stress matrix and the configuration matrix of G(p) respectively. Then SP = 
and Se = 0, where e is the vector of all 1 's. Consequently, rank S <n — r — 1. 

Proof. It follows from ([TJ and (O that the jth row of SP is given by 

su(pY+ t *ik(p k ) T = E <%(y-/) r = o. 

k=\,k^i k:(i,k)eE(G) 

Also, e is obviously in the null space of S. Hence, the result follows. 

□ 



2.1 Dimensional and Universal Rigidity In Terms of Stress 
Matrices 



The following theorem provides a sufficient condition for the dimensional rigidity 
of frameworks. 

Theorem 2 (Alfakih |2 1). Let G{p) be a bar framework on n vertices in W for some 
r < n — 2. If G(p) admits a positive semidefinite stress matrix S of rank n — r—l. 
Then G(p) is dimensionally rigid. 

The proof of Theorem [2] is given in Section [4] It is worth pointing out that the 
converse of Theorem|2]is not true. Consider the following framework [2| G(p) on 5 
vertices in R 2 (see Fig|2]i, where the configuration matrix P is given by 



-3 -5 

1 2 
P= 0-1 

2 
4 



and where the missing edges of G are (1,2) and (3,4). It is clear that G(p) is di- 
mensionally rigid (in fact G(p) is also universally rigid) while G(p) has no positive 
semidefinite stress matrix of rank 2. 

The following result, which provides a sufficient condition for the universal rigid- 
ity of a given framework, is a direct consequence of Theorems Q] and |2] 

Theorem 3 (Connelly ||8j[9), Alfakih [2|). Let G(p) be a bar framework on n ver- 
tices in W,for some r < n — 2. If the following two conditions hold: 

1. G{p) admits a positive semidefinite stress matrix S of rank n — r—l. 
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Fig. 2 A dimensionally rigid framework G(p) in K 2 (in fact G(p) is also universally rigid) that 
does not admit a positive semidefinite stress matrix of rank 2. Note that the points p 2 ,p 4 , and p 5 
are collinear; i.e., G(p) is not in general position. 

2. There does not exist a bar framework G(q) in W that is affinely-equivalent, but 
not congruent, to G(p). 

Then G(p) is universally rigid. 

A configuration p (or a framework G(p)) is said to be generic if all the co- 
ordinates of p , . . . ,p n are algebraically independent over the integers. That is, if 
there does not exist a non-zero polynomial / with integer coefficients such that 
f(p\...,p")=Q. Thus, for a generic framework, Theorem|3]reduces to the follow- 
ing theorem. 

Theorem 4 (Connelly |9|, Alfakih [4|). Let G(p) be a generic bar framework on n 
nodes in W , for some r < n — 2. If G(p) admits a positive semidefinite stress matrix 
S of rank n — r — 1. Then G(p) is universally rigid. 

The proof of TheoremHis given in Section|4] The converse of Theorem|4]is also 
true. 

Theorem 5 (Gortler and Thurston |14|). Let G{p) be a generic bar framework on 
n nodes in W, for some r < n — 2. If G(p) is universally rigid, then there exists a 
positive semidefinite stress matrix S ofG(p) of rank n — r — 1. 

The proof of Theorem [5] given in lfT4l goes beyond the scope of this chapter and 
will not be presented here. 

At this point, one is tempted to ask whether a result similar to Theorem|4]holds if 
the genericity assumption of G(p) is replaced by the weaker assumption of general 
position. A configuration p (or a framework G(p)) in W is said to be in general 
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position if no r + 1 points in p ,p" are affmely dependent. For example, a set of 
points in the plane are in general position if no 3 of them are collinear. The following 
theorem answers this question in the affirmative. 

Theorem 6 (Alfakih and Ye |7|). Let G(p) be a bar framework on n nodes in gen- 
eral position in W,for some r < n — 2. IfG(p) admits a positive semidefinite stress 
matrix S of rank n — r — 1. Then G{p) is universally rigid. 

The proof of Theorem [6] is given in Section [4] The following result shows that 
the converse of Theorem [6] holds for frameworks G{p) where graph G is an (r + 
l)-lateration graph. Such frameworks were shown to be universally rigid in [20|. 
However, it is still an open question whether the converse of Theorem [6] holds for 
frameworks of general graphs. 

A graph G on n vertices is called an (r + \)-lateration graph if there is a permu- 
tation k of the vertices of G, n(l),n(2), . . . , n(n), such that 

• the first (r + 1) vertices, 7f(l), . . . , 7C(r+ 1), induce a clique in G, and 

• each remaining vertex n(j), for j = (r + 2), (r + 3), . . . ,n, is adjacent to (r + 1 ) 
vertices in the set {n(l),7t(2), . . .,7t(j— 1)}. 

Theorem 7 (Alfakih et al [6|). Let G(p) be a bar framework on n nodes in general 
position in W, for some n > r + 2, where G is an (r+ l)-lateration graph. Then 
there exists a positive semidefinite stress matrix S of G(p) of rank n — r — 1. 

The proof of Theorem|7]is given in Section|4] The preceding theorems have been 
stated in terms of stress matrices. The same theorems can be equivalently stated in 
terms of Gale matrices, as will be shown in the next subsection. 



2.2 Dimensional and Universal Rigidity in Terms of Gale Matrices 

Let G(p) be a framework on n vertices in W, r < n — 2, and let P be the configuration 
matrix of G(p). Then the following (r+l) xn matrix 



pT- 




\P 1 . 


.. p n ~ 


e T _ 




1 . 


.. 1 



(4) 



has full row rank since p ,...,p" affmely span W. Note that r <n — 1. Let 

r = the dimension of the null space of i.e., r = n — 1 — r. (5) 



Definition 1. Suppose that the null space of is nontrivial, i.e., r > 1. Any n x r 
matrix Z whose columns form a basis of the null space of is called a Gale matrix 
of configuration p (or framework G(p)). Furthermore, the ith row of Z, considered 
as a vector in W, is called a Gale transform of p' lfl2l . 
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Gale transform plays an important role in the theory of poly topes [ 16 1. It follows 
from Lemma[T]and © that S is a stress matrix of G(p) if and only if 

&>S = 0, and s tj = for all ij (i, j)<^E{G). (6) 

Equivalently, S is a stress matrix of G(p) if and only if there exists an r x r symmetric 
matrix f such that 

S = ZVZ 7 , and s i} = (z'/f Z J = for all ij : i ^ j, (i,j) # E(G), (7) 

where (z') T is the ith row of Z. Therefore, the stress matrix S = Z*PZ T attains its 
maximum rank ofr = «— 1— rif and only if f is nonsingular, i.e., rank W = f, 
since rank S = rank f. 

Then Theorems|2jlll|5]and|6]can be stated in terms of Gale matrices as follows. 

Theorem 8 (Alfakih [2|). Let G(p) be a bar framework on n vertices in W for 
some r < n — 2, and let Z be a Gale matrix ofG(p). If there exists a positive definite 
symmetric matrix f such that 

( z ') T y z j = Ofar each ij : i ? j., (ij) # E(G), 

■ T 

where (z') is the ith row ofZ. Then G(p) is dimensionally rigid. 

Theorem 9 (Connelly (9), Alfakih (4), Gortler and Thurston (TU). Let G{p) be 
a generic bar framework on n nodes in W, for some r < n — 2. Let Z be a Gale 
matrix ofG(p). Then G(p) is universally rigid if and only if there exists a positive 
definite symmetric matrix f such that 

tff¥j = Ofor each ij : i + j, (ij) E(G), 

■ T 

where (z!) is the ith row ofZ. 

Theorem 10 (Alfakih and Ye [7 |). Let G(p) be a bar framework on n nodes in 
general position in W, for some r < n — 2. Let Z be a Gale matrix of G(p). Then 
G(p) is universally rigid if there exists a positive definite symmetric matrix f such 
that 

(zfw z j = Ofor each ij : i + j, (ij) E(G), 

■ T 

where (z l ) is the ith row ofZ. 



3 Preliminaries 

In this section we give the mathematical preliminaries needed for our proofs. In 
particular, we review some basic terminology and results concerning Euclidean dis- 
tance matrices and affine motions. We begin with notation. 
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Throughout this chapter, | . 1 1 denotes the Euclidean norm. |C| denotes the cardi- 
nality of a finite set C. We denote the node set and the edge set of a simple graph 
G by V(G) and E(G) respectively. y„ denotes the space ofnxn real symmetric 
matrices. Positive semi-definiteness (positive definiteness) of a symmetric matrix A 
is denoted by A >z (A >- 0). For a matrix A in y„, diag(A) denotes the n-vector 
formed from the diagonal entries of A. e denotes the vector of all ones in W. A oB 
denotes the Hadamard (element-wise) product of matrices A and B. Finally, the nxn 
identity matrix is denoted by /„; and denotes the zero matrix or the zero vector of 
the appropriate dimension. 



3.1 Euclidean Distance Matrices (EDMs) 

An nxn matrix D = (du) is said to be a Euclidean distance matrix (EDM) if 
and only if there exist points p > ,..., p n in some Euclidean space such that dij = 
\\p' — pi\\ 2 for all i, j = 1, . . . ,n. The dimension of the affine subspace spanned by 
p 1 , . . . ,p" is called the embedding dimension of D. 

It is well known lfTTl[T5l[T8l[T9i that a symmetric n x n matrix D whose diagonal 
entries are all zeros is EDM if and only if D is negative semidefinite on the subspace 

M:={ieR":/i = 0}, 

where e is the vector of all l's. 

Let V be the BX(n-l) matrix whose columns form an orthonormal basis of M; 
that is, V satisfies: 

V T e = 0, V T V =/„_i. (8) 

Then the orthogonal projection on M, denoted by J, is given by J := VV T = I„ — 

ee J jn. 

Recall that S^ n -\ denotes the subspace of symmetric matrices of order n— 1 and 
let y H ^{Ae,y„: diag(A) = 0}. Consider the linear operator & v : y H -> ^ n -\ 
such that 

,? v {D):=- l -V T DV, (9) 
Then we have the following lemma. 

Lemma 2 (|5'|). Let D e ,y H . Then D is a Euclidean distance matrix of embedding 
dimension r if and only if ,%/ (D) y and rank ,%/ (D) = r. 

Let .Jf v : S" n -\ -> y H defined by 

Jf v {X) :=&mg(VXV T )e r + e{&ag(VXV T )) T -2VXV 7 '. (10) 

Then it is not difficult to show that the operators 3y and j£y are mutually inverse 
0. Thus, Lemma|2]implies that D in S?h is an EDM of embedding dimension r if 
and only if D = Jtv(X) for some positive semidefinite matrix X of rank r. 
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Lemma[2]is used in the following subsection to characterize the set of equivalent 
frameworks. 



3.2 Characterizing Equivalent Bar Frameworks 

Since all congruent frameworks have the same EDM, (or equivalently, the same 
projected Gram matrix), in the rest of this chapter we will identify congruent frame- 
works. Accordingly, for a given framework G(p) we assume without loss of gener- 
ality that the centroid of the points p , . . . ,p n coincides with the origin; i.e., P T e — 0, 
where P is the configuration matrix of G(p). 

Let D = (djj) be the EDM generated by framework G(p) in W and let P be 
the configuration matrix of G(p) defined in ||3}. Let X = 3fy(D), or equivalently, 
D = Jffy (X ) ; and let B = PP T be the Gram matrix generated by the points p 1 , . . . , p" . 
Clearly, B is positive semidefinite of rank r. Observe that 

dij= W-AW 

= (p i ) T p i + (p j fp j -2(p i ) T pj, 
= (PP r ),, + (PP r ) ;7 -2(PP r ) ;7 . 



Therefore, 

D = diag(B)e r + e(diag(B)) r - 2B = JfT v (X). 

Hence, 

B = VXV T , andX = V T BV = V T PP T V. (11) 

Furthermore, matrix X is (n — 1) X (n — 1) positive semidefinite of rank r. Accord- 
ingly, X is called the projected Gram matrix of G(p). 

Now let G(q) in W be a framework equivalent to G(p). Let D q and D p be the 
EDMs generated by G(q) and G(p) respectively. Then H o D q = H o D p where H is 
the adjacency matrix of graph G. Thus, 

Ho(D q -D p )=Hoje v (X q -X p )=0, (12) 

where X q and X p are the projected Gram matrices of G(q) and G(p) respectively. 

Let E'i be the n x n symmetric matrix with l's in the ijth and jith entries and 
zeros elsewhere. Further, let 

M ij := & v (E ij ) = -^V T E ij V. (13) 

Then one can easily show that the set {M ,; : i ^ j, ^ E(G)} forms a basis for 
the null space of H o J£y. Hence, it follows from (fT2l that 
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x 9~ x p= E y>r x,u - (W) 

for some scalars yy. Therefore, given a framework G(p) in R'~, the set of projected 
Gram matrices of all frameworks G(q) that are equivalent to G(p) is given by 

{X:X=X p + £ v„.Vf ? " 0}. (15) 

The following lemma establishes the connection between Gale matrices and pro- 
jected Gram matrices. 

Lemma 3 (Alfakih 1 1 1). Let G(p) be a bar framework in W and let P and X be the 

configuration matrix and the projected Gram matrix ofG(p) respectively. Further, 
let U and W be the matrices whose columns form orthonormal bases for the null 
space and the column space ofX. Then 

1. VU is a Gale matrix ofG(p), 

2. VW = PQfor some r x r non-singular matrix Q. 

Proof. It follows from CCD that XU = V T PP T VU = 0. Thus P T VU = 0. Hence, 
VU is a Gale matrix of G(p) since obviously e T VU = 0. 

Now, (VW) T VU = 0. Thus VW = PQ for some matrix Q since P T e = 0. More- 
over, Q is nonsingular since rank PQ = r. 

□ 



3.3 Affine Motions 

Affine motions play an important role in the problem of universal rigidity of bar 
frameworks. An affine motion in W is a map / : W — > W of the form 

f(p i )=Ap i + b, 

for all p' in W, where A is an r x r matrix and b is an r- vector. A rigid motion is an 
affine motion where matrix A is orthogonal. 

Vectors v 1 , . . . , v'" in W are said to lie on a quadratic at infinity if there exists a 
non-zero symmetric r x r matrix <P such that 

(VfcpV = 0, for all i = 1,.. .,m. (16) 

The following lemma establishes the connection between the notion of quadratic 
at infinity and affine motions. 

Lemma 4. ( Connelly MOV ) Let G{p) be a bar framework on n vertices in W. Then 
the following two conditions are equivalent: 
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1. There exists a bar framework G(q) in W that is affinely-equivalent, but not con- 
gruent, to G{p), 

2. The vectors p' — p J for all (i, j) £ E(G) lie on a quadratic at infinity. 

Proof. Suppose that there exists a framework G(q) in W that is affinely- 
equivalent, but not congruent, to G(p); and let q' = Ap' +b for all i = 
Then (q> - q j ) T (q l - q j ) = (p' - p j ) T A T A(p i - pi) = {p< - p j ) T (p l - p j ) for all 
(i,j) £ E(G). Note that matrix A is not orthogonal since G(q) and G(p) are 
not congruent. Therefore, (p' — pi) T <P(p' — p J ) = for all (z',j) £ E(G), where 
<P = I r —A T A is a nonzero symmetric matrix. 

On the other hand, suppose that there exists a non-zero symmetric matrix <P such 
that 0* - P i ) T ^{p i - P j ) = 0, for all (ij) £ E(G). Then l r 8<P >- for suffi- 
ciently small 8. Hence, there exists a matrix A such that /, — 8<P = A T A. Note that 
matrix A is not orthogonal since <P is nonzero. Thus, (p l — p J ) (I r —A T A)(p' — 
p-') = for all (i,j) £ E(G). Therefore, there exists a framework G(q) in W that 
is equivalent to G(q), where q' =Ap' for all i = 1, . . . ,n. Furthermore, G(q) is not 
congruent to G(p) since A is not orthogonal. 

□ 

Note that Condition 2 in Lemma |4] is expressed in terms of the edges of G. An 
equivalent condition in terms of the missing edges of G can also be obtained using 
Gale matrices. To this end, let m be the number of missing edges of graph G and 
let y = (yij) be a vector in M m . Let $(y) be the n x n symmetric matrix whose /j'th 
entry is given by 

&r v \ _ / yij if i + j and (h j)<£E{G), 
g Wh-\ o Otherwise. 

Then we have the following result. 

Lemma 5. (Alfakih [4]) Let G(p) be a bar framework on n vertices in W and let Z 
be any Gale matrix of G(p). Then the following two conditions are equivalent: 

1. The vectors p' — p 1 for all £ E{G) lie on a quadratic at infinity. 

2. There exists a non-zero y = (yij) £ R m such that: 

V T <?(y)Z = 0, (18) 

where V is defined in 0. 

Proof. Let P be the configuration matrix of G(p), and let U and W be the 
matrices whose columns form orthonormal bases for the null space and the column 
space of X, the projected Gram matrix of G(p). Then by Lemma|3]we have 

(p i -p j ) T 0{p i -p J ) = {p'f <V + (p j f <Pp j ~2(p i ) 1 cp p i 
= {P<PP T )ii + [P<PP T )jj - 2{P<PP T )ij 
= (VW0'W T V T )n + {VW0'W T V T )jj - 2(VW0'W T V T )ij 



= jev(w<p'w T ) 
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where <£' = Q<t>Q T for some nonsingular matrix Q, and where J£y is defined in 
CO]). 

Therefore, p' — p J for all G E(G) lie on a quadratic at infinity if and only if 
there exists a nonzero matrix <£' such that H o J£y (W<P'W T ) = 0. But since the set 
{M'-i : i ^ j, ^E(G)} forms abasis for the null space of Ho Jty, it follows that 
vectors p' — pi for all G E(G) lie on a quadratic at infinity if and only if there 
exists a nonzero rxr matrix <P' and a nonzero y = (yij) in R m such that 



(19) 

Next we show that ( fT9b is equivalent to (fT8T l. Suppose there exists a nonzero y that 
satisfies ( fT9] l. Then by multiplying ([19} from the right by U we have that y also 
satisfies ( TT~8T >. Now suppose that there exists a nonzero y that satisfies ( fTST l. Then 



w T 

U T 



'-20' 0" 




~W T ~ 







U T 



= -2W0'W T . 
Thus y also satisfies (TT9b and the result follows. 



□ 



3.4 Miscellaneous Lemmas 

We conclude this section with the following lemmas that will be needed in our 
proofs. We begin with the following well-known Farkas Lemma on the cone of 
positive semidefinite matrices. 

Lemma 6. Let A 1 , . . . ,A k be given n x n symmetric matrices. Then exactly one of 
the following two statements hold: 

1. there exists Y >- such that trace (A'Y) = Ofor all i = 1 , . . . , k. 

2. there exists x = (.*,-) G M. k such that x\A l H Yx k A k >_ 0, ^ 0. 

Proof. Assume that statement 1 does not hold, and let 

if = {Y G y„ : trace (A'Y) = 0, for all i = 1,. . .,k}. 

Then the subspace if is disjoint from the interior of the cone of n x n posi- 
tive semidefinite matrices. By the separation theorem ifTTl page 96], there ex- 
ists a nonzero symmetric matrix such that trace (&Y) — for all Y G if and 
trace (0C) > for all C> 0. Therefore, b and = £f =1 x,^' for some nonzero 
x = (jc/) G R . Hence, statement 2 holds. 
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Now assume that statements 1 and 2 hold and let = x\A l H Vx^A k . Then on 

one hand, trace (&Y) > 0; and on the other hand trace (0Y) = Y%=\ */trace {A'Y) = 
0, a contradiction. Hence, the result follows. 

□ 

The following lemma shows that Gale matrices have a useful property under the 
general position assumption. 

Lemma 7. Let G(p) be a bar framework on n nodes in general position in W and 
let Z be any Gale matrix of G(p). Then any r x f sub-matrix ofZ is nonsingular. 

Proof. Assume r < r. The proof of the case where f > r + 1 is similar. Let Z' 
be any r xr sub-matrix of Z, and without loss of generality, assume that it is the 
sub-matrix defined by the rows f + 1 , r+ 2, . . . , 2f, Then, Z' is singular if and only if 
there exists a nonzero | G 9T such that Z'| = 0. Clearly, Z% is in the null space of 
& '. Furthermore, Z'% = if and only if the components (Z^) f+ i = (Z<^) ?+ 2 = ■ • • = 
(Z| )2p = 0. Now since ZE, ^ 0, this last statement holds if and only if the following 
r + 1 points p l ,p 2 ,... ,p r ,p 2r+l , . are affmely dependent; i.e., G(p) is not in 
general position. 

□ 



4 Proofs 

In this section we present the proofs of the theorems stated in Section|2] 



4.1 Proof of Theorem^ 

Let G(p) be a given framework on n vertices in IR r for some r < n — 2. Clearly, if 
G(p) is universally rigid then G(p) is dimensionally rigid and there does not exist a 
framework G(p) in W that is affmely-equivalent, but not congruent, to G(p). 

To prove the other direction, let X p be the projected Gram matrix of G(p). Let 
Q = [W U] be the orthogonal matrix whose columns are the eigenvectors of X p , 
where the columns of U form an orthonormal basis for the null space of X p . 

Now suppose that G(p) is not universally rigid. Then there exists a framework 
G(q) in M s , that is equivalent, but not congruent, to G(p), for some s: 1 < s < n — 1. 
Therefore, there exists a nonzero y in R™ such that X(y) =X p + ^(y) h where 
■^{y) = Y^(i.i)<£EyijM' ] ■ Now for a sufficiently small positive scalar 8 we haveQ. 

X(ty) = X p + Jt{tf) h 0, and rank(Z(fy)) = rank(X p + Jt(ty)) > r, (20) 

for all t : < t < 8. But, 



the rank function is lower semi-continuous on the set of matrices of order n — 1 
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Q T {X p + jK{ty))Q = 



A+tW T ^(y)W tW T ^{y)U 
tU T .J((y)W tU T JZ{y)U 



to, 



where A is the r x r diagonal matrix consisting of the positive eigenvalues of 
X p . Thus U T ^£(y)U t and the null space of U T JZ{y)U C the null space of 
W T ^(y)U. 

Therefore, if rank (X(to$)) > r + 1 for some < ?o < 5 we have a contradiction 
since G(p) is dimensionally rigid. Hence, rank (X(ty)) = r for all t : < t < 8. 
Thus, both matrices U 1 m($)U and W 1 \j£{y)U must be zero. This implies that 
^{y)U = i.e., V 7 S(y)Z = which is also a contradiction by Lemma|5] There- 
fore, G(p) is universally rigid. 

□ 



4.2 Proof of Theorem^ 

Let G(p) be a given framework on n vertices in W for some r < n — 2 and let Z 
be a Gale matrix of G(p). Let X p be the projected Gram matrix of G(p), and let 
Q = [W U] be the orthogonal matrix whose columns are the eigenvectors of X p , 
where the columns of U form an orthonormal basis for the null space of X p . 

Assume that G(p) admits a positive semidefinite stress matrix S of rank n — r — 1 . 

■ T 

Therefore, there exists a positive definite symmetric matrix f such that (z' ) f z 1 = 
for all ij : i ^ j, £E{G). Hence, by lemma|6] there does not exist y = (yy) € W 

■ T ■ T 

such that Y*ij:i^j,(ij)^E(G)yij(z'(z J ) +Z J (z') ) is a non zero positive semidefinite 
matrix. But z'(z j ) T + z j (z'f = Z T E ij Z. Thus, there does not exist y = (y y ) G M* 
such that Z 7 S(y)Z is a nonzero positive semidefinite matrix. Hence, there does not 
exist y — £ M* such that U T ^(y)U is a nonzero positive semidefinite matrix. 

Now assume that G(p) is not dimensionally rigid then there exists a nonzero y 
such that X = X p + J({y) ^OandrankX > r+ l.But 



e r (x p +^( y ))e: 



U T Jt\y)W U T .£{y)U 



^0. 



Since A +W T ^(y)W is r x r, it follows that U T ^(y)U is a nonzero positive 
semidefinite, a contradiction. 

□ 



4.3 Proof of Theorem^ 



We begin with the following lemma. 
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Lemma 8. (Connelly HlOV ) Let G(p) be a generic bar framework on n vertices in 
W. Assume that each node of G has degree at least r. Then the vectors p' — p 1 for 
all (i,j) E E(G) do not lie on a quadratic at infinity. 

Now let G(p) be a generic bar framework on n vertices in R r . If G(p) admits 
a positive semidefinite stress matrix of rank n — r—1, then each vertex of G has 
degree at least r+l [2 Theorem 3.2]. Thus Theorem[4]follows from Lemmas|4]and 
[8]and Theorem[3] 



4.4 Proof of Theorem® 

The main idea of the proof is to show that Condition 2 of Lemma [5] does not hold 
under the assumptions of the theorem. The choice of the particular Gale matrix to 
be used in equation ( TT~8T > is critical in this regard. The proof presented here is that 
given in J7] . 

Let N{i) denote the set of nodes of graph G that are non-adjacent to node z; i.e., 

N(i) = {j £ V(G) : ./ + i and $ E(G)}, (21) 

Lemma 9. Let G(p) be a bar framework on n nodes in general position in W, r < 
n — 2. Assume that G(p) has a stress matrix S of rank n — 1 — r. Then there exists a 
Gale matrix Z of G(p) such that zij — Ofor all j = 1, . . . ,r and i G N(j + r + 1 ). 

Proof. Let G(p) be in general position in W and assume that it has a stress 
matrix S of rank r=(n—\—r). Let Z be any Gale matrix of G(p), then it follows 
from (01 that S = Z*FZ T for some non-singular symmetric r x r matrix f . Let us 
write Z as: 

"Zi' 



z = 



(22) 



Z 2 

where Z2 is f x r. Then it follows from Lemma|7]that Z2 is non-singular. Now let 

Z=(zij)=ZVZ2 T . (23) 
Then Z is a Gale matrix of G(p) since both f and Z2 are non-singular. Furthermore, 

S = Z¥Z T = Zf [z\ Zj] = [ZWZj Z\. 

In other words, Z consists of the last f columns of S. Thus Zij — Sij+ r +i. It fol- 
lows by the definition of S that su+ r +i — for all i, j such that i ^ (j + r+l) and 
(i,j + r+ 1) $E(G). Therefore, z u = for all j = 1 , . . . , f and i G N(j + r + 1 ). 

□ 

Lemma 10. Let the Gale matrix in ( 17^1 ) be Z as defined in rf23l ). Then the system of 
equations M&\ is equivalent to the system of equations 



S{y)Z = 0. 



(24) 
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Proof. System of equations ( fTSI ) is equivalent to the following system of equa- 
tions in the unknowns, yy (i ^ j and £ E(G)) and | = (^) e M r : 

<g(y)Z = et; T . (25) 

Now for /' = 1 , . . . , r, we have that the (j + r + 1 , j)th entry of (y)Z is equal to ^j. 
But using (TTTb and Lemma [9] we have 

n 

i=l iuG/VO'+r+l) 

Thus, t, — and the result follows. 

□ 

Lemma 11. Let G{p) be a bar framework on n nodes in general position in W, 
r < n — 2. Assume that G(p) has a positive semidefinite stress matrix S of rank 
r = n — 1 — r. Then there does not exist a framework G(q) in W that is affinely- 
equivalent, but not congruent, to G(p). 

Proof. Under the assumption of the lemma, we have that deg(z") > r + 1 for all 
i G V(G), i.e., every node of G is adjacent to at least r+1 nodes (for a proof see (2] 
Theorem 3.2]). Thus 

\N(i)\ <n-r-2 = r-lforalH GV(G). (26) 

Furthermore, it follows from Lemmas l9l [TOl and 151 that the vectors p' — p ; for all 
£ E(G) lie on a quadratic at infinity if and only if system of equations d24l i has 
a non-zero solution y. But (|24| > can be written as 

Y_ yijZ J = 0, for/= l,...,n, 

where (z') r is the ith row of Z. Now it follows from ( f26b that yy = for all (i,j) £ 
E(G) since by Lemma|7]any subset of {z l ,...,£"} of cardinality < r — 1 is linearly 
independent. 

Thus system (l24l i does not have a nonzero solution y. Hence the vectors p' — p>, 
for all (i,j) 6 E(G), do not lie on a quadratic at infinity. Therefore, by LemmaH] 
there does not exist a framework G(q) in W that is affinely-equivalent, but not con- 
gruent, to G{p). 

□ 

Thus, Theorem|6]follows from LemmafTTIand Theorem|3] 
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4.5 Proof of Theorem^ 

The proof of Theorem [7] is constructive, i.e., an algorithm is presented to construct 
the desired stress matrix. The proof presented here is a slight modification of that 
given in (6). 

Let G(p) be a framework on n vertices in general position in W, n > r + 2, and 
let Z be a Gale matrix of G(p). An n x n symmetric matrix S that satisfies 

— 0, or equivalently S = Z*¥Z T for some symmetric matrix f, 

is called a pre-stress matrix, where & is defined in (0). Thus, it follows from © 
and (0 that S is a stress matrix of G(p) if and only if S is a pre-stress matrix and 
Sij = for all ij : i ^ j, (i,j) E(G). 

Clearly, S" — ZZ T is a positive semidefinite pre-stress matrix of rank r = n — 
r — 1, If S" satisfies s"j — for all ij : i ^ j, $ E(G), then we are done since 
S" is the desired stress matrix. Otherwise, if S" is not a stress matrix, we need to 
zero out the entries which should be zero but are not, i.e., the entries s"- ^ 0, i ^ j 
and £ E(G). We do this in reverse order by column (row); first, we zero out 
the entries s" n ^ 0, for i < n and (i,n) $ E(G), and then do the same for columns 
(rows) (n — 1), (n — 2), . . . , (r + 3). This "purification" process will keep the pre- 
stress matrix positive semidefinite and maintain rank n — r — 1 . 

Let G be an (r + l)-lateration graph with lateration order 1,2, ...,n; i.e., the 
vertices, 1,2, ...,r + 2, induce a clique in G, and each remaining vertex k, for 
k = r + 3, . . . ,n, is adjacent to (r + 1) vertices in the set {1,2, 1}. Let 



We first show how to purify the last column (or row) of S" — ZZ T . Let Z" denote 
the sub-matrix of Z obtained by keeping only rows with indices in N'(n) U {«}. Then 
Z™ is a square matrix of order f = n — r — 1. Furthermore, by Lemma [7] it follows 
that Z" is nonsingular. Let b" denote the vector in W such that 



N'(k) = {ieV(G):i<k and (i,k) &E(G)}. 



(27) 



Then for k = r + 3 , . . . , n, 



\N'(k)\=k-r-2, 



(28) 




1 if i = n. 




Now let |„ G M' be the unique solution of the system of equations 



Lemma 12. Let S' 



71—1 



S" +Z^ n T Z T = Z(I + ^,, T )Z T . Then 



1. S" 1 is a pre-stress matrix ofG(p), i.e., &S n 1 = 0. 

2. S"~ l y and the rank ofS' 1 ^ 1 remains n — r — 1. 

3. s'l- x = 0foralli:i <n, (i,n) &E(G). 
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Proof. The first statement is obvious. The second statement follows since / + 
t,n<^n > 0. The third statement is also true by construction. For all i < n, (i,n) 
E(G), i.e., for all i e N'{n), we have s^ 1 = s'l+b^ = sf n - s n in = 0. 

□ 

We continue this purification process for columns (n — 1), (r + 3). Be- 

fore the feth purification step, we have S >r 0, ,9- 'S k — 0, rank S = n — r — 1, and 

4 = 0, for all ij : i ? j, (ij) # E(G), and for all j = k+ l,...,n. 

Let Z k denote the sub-matrix of Z obtained by keeping only rows with indices in 
N'(k) U {k, k + 1 , . . . ,n}. Then Z k is a square matrix of order r = n — r — 1 . Further- 
more, by Lemma|7] it follows that Z k is nonsingular. Let b k denote the vector in W 
such that 

<-s\ k iU&N'{k), 
b k =l 1 if i = k, 

[ if i = k+\,...,n. 

Now let §j G M r be the unique solution of the system of equations 

Z% = b k . 

The following lemma shows results analogous to those in Lemma [T2l for the 
remaining columns. 

Lemma 13. Lets' 1 - 1 = S k + Z^Z T . Then 

1. S k ~ l is a pre-stress matrix of G(p), i.e., £5 g S k ~ 1 = 0. 

2. S k ~ l >z and the rank ofS k ~ l remains n — r — 1 . 

3. s k J l —Q for all i : i < j, (i,j)(£E(G) andforallj = k,...,n. 

Proof. The proof of the first two statements is identical to that in Lemma [T2] The 
third statement is again true by construction. For each i < k, (i,k) ^ E(G), i.e., for 
all i 6 N'(k), we have 

s ik = s ik+ b i b k = s ik- s ik = Q - 

Furthermore, for j = k + 1 , . . . ,n, the jth column (or row) of Z ^^Z T has all zero 
entries, which means that the entries in the y'th column (or row) of S k ~ l remain un- 
changed from S k . 

□ 

Proof of Theorem 

The matrix 

S r+2 = S'-+ 3 +Zt; r+ & + 3Z T = Z(I + &g + ■ ■ ■ + U^J +i )2\ 

obtained at the "(r + 3)th" step of the above process, is by Lemmas fl~2l and [T3l a 
positive semidefinite pre-stress matrix of rank n — r — 1 . Furthermore, s'^ 2 = for 
all ij : i ^ j, (i, j) E (G) and for all /' = r + 3, r + 4, . . . , n. But since the vertices 
1 , 2, . . . , r + 2 induce a clique in G, it follows that 
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sJ+ 2 = for all ij-.i^j, (iJ)^E(G). 

Hence, S — S r+1 is a positive semidefinite stress matrix of G(p) of rank n — r — 1; 
i.e., S r+2 is the desired stress matrix. 
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